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Abstract 

We explicitly show that split-universal extra dimension (split-UED), a recently suggested exten- 
sion of universal extra dimension (UED) model, can nicely explain recent anomalies in cosmic-ray 
positrons and electrons observed by PAMELA and ATIC/PPB-BETS. Kaluza-Klein (KK) dark 
matters mainly annihilate into leptons because the hadronic branching fraction is highly suppressed 
by large KK quark masses and the antiproton flux agrees very well with the observation where no 
excess is found . The flux of cosmic gamma-rays from pion decay is also highly suppressed and 
hardly detected in low energy region < 20 GeV). Collider signatures of colored KK particles 
at the LHC, especially q\q\ production, are studied in detail. Due to the large split in masses of 
KK quarks and other particles, hard pt jets and missing Et are generated, which make it possible 
to suppress the standard model background and discover the signals. 
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I. INTRODUCTION 



Recent observations by PAMELA Q and ATIC/PPB-BETS |2|,|3| consistently suggest a 
new primary source of energetic cosmic electrons and positrons in the neighborhood of our 
solar system, since high energy electrons/positrons lose their energy quickly within about 1 
kpc. These resutls have attracted enormous attentions and stimulated many interpretations 
of the primary electric source from both astrophysics and particle physics, for example, 
the nearby pulsars 4] and dark matter decay 5] or annihilation The existing data of 
cosmic electrons and positrons cannot distinguish among these possibilities, and we expect 
that future data from cosmic-ray and accelerator experiments such as Fermi and Large 
Hadron Collider (LHC) will give us some clues on which interpretation is more promising 
and ultimately correct. 

In this paper, we consider one of the most attractive models of dark matter, universal 
extra dimension (UED) with its minimal realization (mUED)^, 8] and a recently suggested 
variety split-UED 9|. When the reflection symmetry (dubbed Kaluza-Klein parity) about 
the mid point of the orbifold extra dimension is exact, the lightest Kaluza-Klein particle 
(LKP), which is odd under the parity operation, is absolutely stable and is a good candidate 
of dark matter [ji]] . It turns out that in most cases of universal extra dimension scenarios the 
LKP is the first Kaluza-Klein (KK) photon (i?i), and a pair of charged leptons is produced 
in the annihilation process of B\. Interestingly enough, the mass of dark matter suggested 
by the ATIC/PPB-BETS data (~ 600 — 700 GeV) exactly coincides with the one that gives 
the right relic density for the LKP dark matter after taking co-annihilation channels into 



account 



ill 



On the other hand, the measured data of the cosmic ray antiproton-to-proton flux ratio 
between 1 and 100 GeV by PAMELA follow the expectations from secondary production 
and strongly constrain contributions from dark matter particle annihilation 12J. Therefore 
it is required to naturally reduce antiproton flux in any dark matter model. Split-UED is 
promising in the sense that production of antiproton is quite suppressed. The cross section 
for dark matter annihilates into a quark pair is 

/ \ 2 // 2 i 2 \2 

[a Bl B^qqv) oc m B J(m Bl + m qi ) 

where m Bl and m qi are the masses of £>i and the first KK quark (qi), respectively, and it 
can be significantly suppressed by increasing the masses of KK quarks, i.e. splitting the 



KK quarks from other particles. However we should keep in mind that predictions of the 
observed antiprotons from dark matter annihilation can vary quite a lot for different adopted 
diffusion models as we will see later. As hadronic production is suppressed in split-UED, we 
would also expect much less production of cosmic gamma-ray whose main source is decay of 



pions. This prediction will be tested by forthcoming data from the Fermi experiment 13]. 



The split spectrum of Kaluza-Klein particles is realized by introducing a bulk mass term 
for quarks in split-UED keeping the KK parity exact 9| 1 . One immediate consequence 
is that the 5D wave functions of quark fields are quasi-localized at the boundaries and 
induce violation of the translational symmetry along the extra dimension. Accordingly the 
KK number conservation is violated in the quark sector, and KK-even gauge bosons could 
directly couple to the zero mode quarks at the tree level, which is forbidden in mUED. 
KK-even bosons could be copiously produced at colliders and it is promising to discover 
them, for example, by the 2nd KK gluon resonance at the LHC 9j. Another consequence 
is the heavy KK quark decay to the SM quark and KK boson, which will generate harder 
jets with large missing momentum compared to that in mUED due to a large mass splitting 
between KK quarks and other KK particles. Such differences may give us some handle to 
tell split-UED apart from mUED at the LHC. 

This paper is organized as follows. In Section [Til we review split-UED and show how 
the spectrum and the interactions are modified compared with those in mUED. Note that 
the renormalization group (RG) running is taken into account in our spectrum calculation. 
We then calculate the annihilation cross section to leptons and quarks with varying "bulk 
mass parameter", /i, which is the only new extra parameter in split-UED. In Section IITT1 
we estimate the signals of cosmic-ray positrons, electrons, and antiprotons from LKP dark 



matter annihilation in both mUED 15j and split-UED, and compare our predictions with 
recent experimental data. Moreover, the diffusive gamma-ray is studied as well. In Section 
\IV\ we consider the productions of KK quarks and gluons at the LHC, and we also discuss 
the possibility of discovering these signals and mass reconstructions. We finally conclude 
our studies in Section 



1 In warped geometry KK parity could be kept in a different way [14j . 
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II. MASSES AND COUPLINGS IN SPLIT-UED 



Our set-up is based on the all SM fields are "universally" propagating in one extra dimen- 
sion, which is compactified on an orbifold, S 1 /Z 2 , with the boundary points y = ±L (the 
orbifold radius R = 2L/tt). The extension we have introduced is an scalar background 
that couples to each colored Dirac fermions ^/ q (x IJ ', y) in the 5D bulk where the corresponding 
standard model quark field q (=SU(2) singlets u c R , d c R and doublet Ql) resides on the chiral 
zero mode after orbifold projection. A common Yukawa coupling A(= Xq l = \ u c R = Ad c ) 
is assumed for all quarks so that any further source of flavor violation beyond the CKM 
matrix is automatically avoided and the model becomes more predictive 2 . Here we neglect 
the family indices but one can easily extend the model to include more generations. Finallly 
the 5D action is given by the form 

^spiit-UED = J d 4 x J dy £ mU ED - \<&(y)^q(x, y)ty g (x, y) , (1) 

where the summation runs for all quark fields q, and £ m uED contains the usual mUED terms. 

If we choose an odd background profile $(— y) = —&{y) under the inversion symmetry 
about the middle point y = 0, we can see that the KK parity, which transforms the fields as 
<3>(x, y) — > — y), ^ q (x, y) — »• ±75 — y), is still a good symmetry of the Lagrangian. 
For quantitative concreteness we consider the simplest mass profile which respects the re- 
quirement: A($(y)) = /ie(y), where /1 is the bulk mass parameter and e(y) is a step function 
defined as +1 for < y < L and —1 for — L < y < 0. As (A — > 0) or equivalently (/1 — > 0), 
split-UED is continuously reduced to mUED thus we call this limit as "mUED limit" below. 
For gauge bosons, scalars, and leptons (7), their KK decompositions are the same as those 
in mUED, while for quarks, their odd masses will affect the KK decompositions. The de- 
tails of the all KK decompositions in Eq. (CD) are presented in Appendix A. Notice that in 
the case of KK decomposition, we also label the KK parity in addition to the KK number 
n = {n~(n odd), n + (n even)}. 

The tree level mass for the n-th KK fermion is given by 



tree 



^ ^W+K, (2) 



< ee = n/R , (3) 



2 In general the Yukawa coupling can be assigned for any fermion field in the bulk and they can be different 
from each other. 
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if we neglect their masses coming from electroweak symmetry breaking. Here k n is deter- 
mined by the boundary conditions which are given by: 

k n - = — | \i | tan k n L, (4) 
k n + = n/R . (5) 

We choose the minimal universal boundary conditions at the energy cutoff scale A that all 
boundary kinetic terms vanish 3 , and t he g eneral expression for the one-loop correction to 



the 77,-th KK fermion masse is given by 



17| 



tree 



l + ^l^^iV^jlogAi? 



(6) 



G 

where «g = g 2 G j^ and we neglect the Yukawa corrections in this formula 4 . The sum 
is over all SM gauge groups under which the fermion / is charged, and Cf{N) is the 
quadratic Casimir operator of the fundamental representation N in gauge group G, e.g., 
C2 U{N) (N) = (N 2 - 1)/2N and we take Cf (1)y (N) = Yf with Y f being the hypercharge of 
/. If we choose (a u a 2 ,a 3 ) ~ (0.010,0.033,0.094) at Q 2 = (620GeV) 2 , colored particles will 
get bigger corrections (~ 14% for u R , d l R and 16% for Q l L ) than leptons (~ 1.9% for l R and 
3.1% for 11). 

In Fig. [I] we plot the spectra of 1st KK quarks and leptons including 1-loop corrections. 
When fi = 0, i.e. the mUED limit, quarks are quite degenerate with leptons even we take 
the QCD corrections to the masses into account. However when we increase the bulk mass 
parameter fi, the KK quark masses become larger and larger so that we get the split spectra 
as we request. Note that small deviations in quark spectra come from different SU(2) and 
£7(1) charge of u R , d R and Ql. In our study, these differences can be neglected. 

The couplings between KK gauge bosons and fermions determine the most interesting 
phenomenological features of the model. All 4D effective couplings could be calculated by 
integrating out the 5D profiles of the relevant fields along the extra dimension. Since 5D 
profiles are the same for different KK gauge bosons (quarks) with the same KK number, the 
ratio between a general gauge boson- quark- quark coupling gA m -q„-q[ over the SM one is the 



3 One can also consider the non- universal boundary conditions as an extension. See for example, Ref. 16 1. 

4 For the KK top quark, if we consider the large top Yukawa coupling, positive contribution to the tree level 
mass term will cancel the negative one from the radiative corrections, and the overall mass corrections 
will not be affected too much. 
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FIG. 1: (upper plot) The spectra of 1st KK quarks and leptons in the unit 1/R with the given 
bulk mass parameter /i. The lepton masses are independent of fi in our scenario. Here we included 
1-loop correction for all the masses. The most important correction is QCD correction and it 
is roughly 14 ~ 16% for quarks. For leptons the correction is small (~ 2%). (lower plot) The 
magnified figure of the upper one in fiR = (3.60,3.90). 577(2) doublets (tL,Qi) get the larger 
1-loop correction than singlet quarks (ur,cIr) get. 

same for fixed KK numbers, where m, n and / are KK number of A, q and q', respectively. 
We introduce a function Q mn i{^R) to parameterize how such a ratio depends on the bulk 
mass parameter /x times the radius R: 

9A m -g n - qi /90 = Gmnli^R) , (7) 

where g is the SM gauge interaction between A , q and q' Q . Among those couplings, the 

6 



0.4 







2 3 



FIG. 2: The coupling constant ratio gA 1 ~q 1 ~q a /9o = Giw vesus fiR. The dependence comes from 
the change of the 5D profiles. Notice that the coupling ratio first increases to its maximal value 
around 1.14 for \xR ~ 1 and then decreases monotonically. 

most interesting one in phenomenology is the coupling between the 1st KK gauge boson 
Ai, the 1st KK quark qi, and SM quark q , which is crucial on LKP pair annihilation and 
productions of q\ and the 1st KK gluon g\ at the LHC. The dimensionless function Giioi^R) 
is calculated from the 5D profiles of Ai, qi and q given by Eq. (1571) . (1551) . (1551) . and (l4"Tj) in 
Appendix A, and is plotted in Fig. [2] as a function of fiR. At fi = limit we get the mUED 
result, which is the same as the SM coupling. As \x gets larger, the 5D profiles are localized 
so that the resultant overlap among the fields get smaller. 

Taking care of all the corrections and effects of /i we can now calculate the branching 
fractions that LKP pair annihilates into quarks and leptons. At the non-relativistic limit 
{(3 = y/l — 4mg i / s — > 0), the cross section for LKP pair annihilates into fermions can be 
well approximated as: 



W v )b 1 b 1 - 



§Txm 2 Bi 
N c (Y f \ H 



+ 0{v) 



(8) 



Yf ) is introduced with the color factor 



where a useful numerical term Cf 
3(1) for quarks(leptons). The mass ratio between the LKP and KK fermion is parameterized 
by a parameter rj defined as Tf = rrif l /mB 1 and we ignore the small mass corrections coming 




FIG. 3: The ratio of LKP pair annihilation into quarks and charged leptons vesus [iR. The mUED 
limit corresponds to [i = and the ratio of cross section to quarks to leptons is about 47% but it 
quickly goes down below 10% as fi > 1.3/R. The shaded region (a qq jall > 10%) is not preferred 
by PAMELA antiproton data. 

from SU(2)l, U(1)y gauge group so that rrif 1 = mf L1 ~ m f R1 - In mUED limit, the ratio 
between productions of the lepton and the quark is mainly controlled by the ratio of Cf from 
the hypercharges and color factor of fermions. For leptons and quarks the corresponding 
terms are ^ C\ ~ 3.18 and ^ q C q ~ 1.90, respectively. Obviously the mUED does have a 
sizable hadronic annihilation channel and a carefully study of antiproton flux is important. 
In the case of split-UED, one can easily adjust r q to suppress the hadronic annihilation cross 
section. In Fig. [3] we plot the ratio of annihilation cross section into quarks and charged 
leptons cr qq /aii with /i_R in the range of (0, 5). We find the ratio goes down very quickly as 
the bulk mass parameter goes larger since a ~ 1/(1 + f 2 ) 2 . By turning on the bulk mass 
parameter we can easily avoid the antiproton excess in the cosmic ray detection experiments. 
We also note that the soft gamma-ray from decay of hadrons, especially from pions, is under 
control in this set-up. 
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TABLE I: The branching fraction of the LKP pair annihilates into various final states, for different 
q± mass, and different bulk mass parameter fj,. The numbers shown are summed over generations. 
We neglect the TV-boson final state because the SU(2) component of LKP is too small. Here we 
fix the LKP mass to be 620 GeV and include the radiative corrections to all the masses according 



to Ref. U, 



H (GeV) 





200 


400 


600 


800 


1000 


M qi (GeV) 


713 


863 


1026 


1198 


1378 


1566 


BR( J B 1 S 1 qq) 


29.4% 


26.4% 


20.6% 


14.3% 


8.9% 


5.2% 


BR( J BiBi -» 11) 


64.3% 


67.1% 


72.3% 


78.2% 


83.0% 


86.5% 


BR( J BiSi -» uu) 


3.8% 


3.9% 


4.3% 


4.6 % 


4.9% 


5.1% 


BR(B 1 B 1 -» tfxj>*) 


2.3% 


2.4% 


2.6% 


2.8% 


3.0% 


3.1% 



III. COSMIC POSITRON, ANTIPROTON AND PHOTON 

LKP dark matter B\ will mainly annihilate into a fermion pair while the W-boson and 
Higgs boson modes are suppressed, as shown in Table [H For the quark-pair final state, due to 
the QCD hadronization process, a bunch of hadrons will be produced and sequentially decay 
into positrons/electrons, protons/antiprotons, photons and neutrinos. Moreover, positrons 
and electrons also originate in the leptonic channels: muon and tau (also generate photons) 
decay and direct production in B\B\ — > e + e~ process. In our calculation, we use PYTHIA 
[bsl for the simulation of the QCD hadronization and the decays of particles. Since these 
stable particles from the dark matter annihilation in our Galaxy will propagate to our solar 
system and be observed as cosmic-ray signals, we estimate the predictions from our model, 
the split-UED, and compared with experimental data in this section. Due to the fact that 
excesses were observed in the positron fraction in PAMELA ]J and the flux of electron plus 
positron in ATIC/PPB-BETS {2, 3] while the ratio of antiproton to proton, p/p, is consistent 
with the astrophysical background, we first explain positron and electron data from the dark 
matter annihilation. Then base on the parameters set by fitting the ATIC/PPB-BETS and 
PAMELA data, we calculate the p/p and photon flux. Since the W-boson and the Higgs 
boson productions are highly suppressed in UED models we are considering, we will not 
take into account their contributions in the following calculations. 
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Model 


5 


A' (kpc 2 /Myr) 


L(kpc) 


M2 


0.55 


0.00595 


1 


MED 


0.70 


0.0112 


4 


Ml 


0.46 


0.0765 


15 



TABLE II: The parameters in the diffusion models which are comparible with B/C and generate 
minimum (M2), median (MED) and maximal (Ml) positron flux. 

A. Positron and Electron 



In this section, we calculate the cosmic-ray positrons and electrons from the -Bi-Bi annihi- 
lation which account for the PAMELA and ATIC/PPB-BETS data. Since the dark matter 
annihilation generates the same amount of electrons and positrons, we only show the formula 
for positron flux below. After being produced from the dark matter annihilation process, the 
positrons will propagate in the magnetic field of the Milky Way. Because the magnetic fields 
are tangled, the motion of the positron can be described by a diffusion equation. Neglecting 
the convection and annihilation in the disk, the steady state solution must satisfy 



K{E) v 2 f e+ (E,r) + ^ [b(E)f e+ (E,r)} + Q(£,f) = 0, 



(9) 



where f e + is the number density of e + per unit kinetic energy, K(E) = Kq(E /GeV) 5 is the 
diffusion coefficient, b(E) = 10~ 16 (-E/GeV) 2 is the rate of energy loss and Q(E,f) is the 
source of producing e + . In the £>i dark matter annihilation case, 



Q(E,r) 



1 /p(f) 

2 \m B , 



fdN(E) ( 



dE 



(10) 



where dN e + / dE is the energy spectrum of e + obtained by using PYTHIA, the index i runs 
over all quark and charged lepton pairs and p(r) is the dark matter profile. In our numerical 
calculations, we adopt an overall boost factor B F and the isothermal halo model [3] which 



is given as 



Phalo(^) 



Po 



(11) 



1 + (r/r c ) 2 

where r = \r\ is the distance from our Galactic center, r c = 3.5 kpc and po is the parameter 
that is adjusted to yield a dark matter local halo density of 0.3GeV/cm 3 [2] in our solar 



Effectively we are using p 2 {r) = Bpp 2 halo (r) 
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system. Considering the diffusion zone of a cylinder with radius R and half-height L, the 
solution for the positron flux at the Earth can be written in a useful form {^(J 



®T(E) = ^f e+ (E,r Q ), (12) 



where c is the speed of light, 



M E,r a ) = f"" 1 ^dE> £ Q„, m (-B')Jo (^) =in (^) exp + (|f ) X 

n,m=l ^ ' L ^ ' 

(13) 

where J^'-* is the energy in a unit of GeV, Jo is the zeroth-order of the first kind Bessel 
function, r ~ 8.5 kpc is the distance from Milky Way center to the Sun, ( n is the n-th zero 
of the function Jo, 

X = g °* 10 " [E 5 ' 1 - E' 5 " 1 ] (14) 

and 

Q "- (E) -M^mSy"^ L dzA " - z) ) Q{E ' r ~>- (15) 



The Ji in the Eq. (115j) is the first-order of the Bessel function and the parameters 5, Kq are 
set in the simulations of diffusion models to coincide with the observed cosmic-ray data, 



especially the Boron to Carbon ratio (B/C) 2lfl . The values for different diffusion models 
we adopted are listed in table HT1 22 ] . 

In addition to e ± fluxes from dark matter annihilation, there exist secondary fluxes 
from interactions between cosmic-rays and nuclei in the interstellar medium. We use the 



approximations of the e and e + background fluxes 23l. |24| 



0.16£ 



-l.i 



7E - 7 

$f- c (£) = — — — GeV-'cm-^ec-^r- 1 , 

e v ; 1 + 110E 1 - 5 + 600E 2 - 9 + 580E 4 - 2 
4 ^770.7 

*S°W = TTgTT500^ GeV-c m -W Sr -, (16) 

where E is in units of GeV. Therefore, the fraction of positron flux and total flux of positron 
plus electron are 

e + + e+ (17) 
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FIG. 4: The flux of electron plus positron from LKP dark matter annihilation using MED diffusion 

nnnn 

model, compared with experimental data [21. la. 126. 1271] . Note that only lepton mode is inclulded in 
the split-UED, since all the KK quarks are decoupled. 



and 



DM 



prim 



(18) 



respectively, where A; is a free parameter which is used to fit the data when no primary 
source of e + flux exists 24j, |25|. For our numerical simulation, we take the mass of dark 
matter, to be 620 GeV. In FigHJ we show the total flux of electron and positron with 



the background from Eq. (1161) and recent experimental data 



as 



271 ]. In order to have 



a better fit for the data, a boost factor of 200 is needed, which is consistent with the value 



(200) chosen in Ref.(^|. Boost factor is known to have origins such as local clumps in dark 



matter profile |28l | (se e a 



force 



30 



31 



so 



291]). Sommerfeld enhancement effect by a long range attractive 



L3J, [33|, EL the Breit-Wigner type resonance effect |35| and possibly many 
more origins. In our case, without assuming a new attractive force or further tuning of mass 
spectra to get the large resonance, we tend to assume that boost factor mainly comes from 
clumps. We also show in FigJH the contributions of quark mode and leptonic mode from 
the dark matter annihilation. The e from quark mode are much softer than that in the 
leptonic mode. This is because that the main source of in quark mode is the hadron 
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comparison of diffusion models 
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FIG. 5: The same as FigHJbut only for split-UED, and different diffusion models, Ml, MED and 
M2, are shown together. 

cascade decay while in leptonic mode there exists a direct production of -Bi-Bi — > e + e~, and 
from ^ and decay are harder. We can easily see that the peak is mainly contributed 
by the leptonic mode and sharp drop-off is due to the direct production of e ± , therefore, 
the predictions of mUED and split-UED are quite similar to each other in the peak region 
since the difference between these two models is just in quark sector. We also compare 
different diffusion models in Figj5]for split-UED case. MED and Ml models produce more 
soft electrons and positrons than the M2 model does, thus the distributions of the former 
two models are flater. For energetic positrons and electrons, all of three models are quite 
the same with each other. 

For positron fraction, we show in Figj6]the predictions of mUED and split-UED, using 
MED diffusion model, with PAMELA data l|. Since there are more soft e + in mUED from 
the quark modes, the distribution is falter compared to that of split-UED. If we adopt M2 
and Ml models, as shown in FigJTJ the fittings seem worse compared the MED case for split- 
UED especially in the energy region around 10 GeV. However, we have decoupled all of the 
higher KK quarks in this plot, i.e. hadronic final states are turned off in the dark matter 
annihilation. With finite mass of a first KK quark, more soft positrons will be produced, 



13 



1 



- m D = 620 GeV, boost factor = 200 

- propagation model: MED 




split-UED ( KK quarks decoupled ) 
mUED 



0.01 - 



10 



100 



1000 



E (GeV) 



FIG. 6: The positron flux fraction predicted from the LKP dark matter annihilation using MED 
diffusion model, compared with PAMELA data 

therefore, the M2 model can also be consistent with the data. For Ml model, we expect 
that by tuning the parameters, e.g. boost factor or the normalization of the background, it 
would be possible to explain the PAMELA data as well. So, we conclude here that the dark 
matter in both mUED and spilt-UED can be the equivalently good source of the excesses 
in PAMELA and ATIC/PPB-BETS data. 

B. Antiproton 

The propagation of antiprotons through the Galaxy is similar to that of positrons. The 
diffusion equation of antiprotons can be described as 



K p v 2 fp(T,r) - — [{sign of z)V c (z) f p (T , r}] - 2hS(z)T ann f p (T,r) + Q(T,r) = 0, (19) 



where fp(T, r) is the number density of antiproton per unit energy, T is the kinetic energy 
of antiproton and K p = K 0p j3(p/GeV) s is the diffusion parameter. V c (z) in the second term 
of Eq. ffl~9]) is related to the convective wind that tends to push antiprotons away from the 
Galactic plane, and is assumed to be a constant . Again, like the case in positron, the values 
of K 0p , S and V c for different diffusion models are set to agree with the observed cosmic-ray 
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FIG. 7: Same as Figj6]but only for split-UED, and different diffusion models, Ml, MED and M2, 
are shown together. 



data, especailly the B/C, and their values are listed in table ITO1 221 ] . The third term of 
Eq. (jT9l) represents the annihilation of p with the interstellar proton in the Galactic plane, 
h is the half-height of plane which is set to be 0.1 kpc in our calculations. The antiproton 



annihilation rate T ann is given as 36] 



Tann = (ji H + 4 2 / 3 n He )a;; n (T)v p , (20) 

where n# ~ 1/cm 3 and % e ~ 0.07n# are the number densities of hydrogen and helium, 
respectively; cr|p n (T) is the annihilation cross section and Vp is the velocity of antiproton. 
The form of &pp n {T) is given by 



661(1 + 0.0115T 
36T -a5 mbarn 



-0.774 



0.948T ' 0151 ) mbarn , for T < 15.5 GeV; 

,forT > 15.5 GeV. 



The solution of the interstellar flux of antiproton in the vicinity of solar system is 



$ IS (T) = gfp(T,rl) 



4rr 2m 



i 2 ^ 



< <7v >i G(T) 



dNp 



[21) 



(22) 
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Model 


5 


K 0p (kpc 2 /Myr) 


L(kpc) 


K(km/s) 


MIN 


0.85 


0.0016 


1 


13.5 


MED 


0.70 


0.0112 


4 


12 


MAX 


0.46 


0.0765 


15 


5 



TABLE III: The parameters in the diffusion models which are comparible with B/C and generate 
minimum (MIN), median (MED) and maximal (MAX) antiproton flux. 

where the dNp/dT is the spectrum of antiproton which we simulate with PYTHIA, and the 
index i runs over all quark final states in the B\ dark matter annihilation processes, and 

Vn(L) 



n,m=l 



V C L 
2K p (T) 



A n smh(S n L/2) 



(23) 



with 

Vn(L) 
and 



drrJ (( n r/R) / dzexp 
J 1 (UJ-K Jo Jo 



V C (L - z) 
2K P (T) 



sinh(S n (L - z)/2)p{r)\ (24) 



A n = 2hT ann + V c + K p (T)S n coth(S n L/2), 



Sn 



KliT) ' R 2 



(25) 
(26) 



However, we have to take solar modulation into account to estimate the flux of antiproton 

"he flux of 

the antiprotons at the top of Earth's atmosphere can, therefore, be written as 

2m p T TOA + T| OA IS 
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1C 


39 


L 



TOAJ 



is, 



(27) 



2m p T ls + T/ s 

where Ttoa = ^is — <Pf with (f)p being the solar modular parameter which we take to be 500 
MV in our calculation. In order to compare with the observed data of antiproton to proton 
ratio, we need to include the astrophysical background. Here we adopt simple fittings for 



background antiprotons $p 9 and protons provided in Ref. [41] 
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(28) 
(29) 



The predicted ratios of antiproton to proton, p/p, for mUED and split-UED are shown in 
Fig. [8] with experimental data 12, 42, 43], 44, 45]. Here we take a common boost factor 



{Bp = 200) which we have obtained from ATIC/PPB-BETS fit. Since the hard electronic 
signals are mainly from local sources within a few kpcs but antiprotons can come from 
distant sources, one may take different values of boost factor for positron and antiproton in 
one's calculation. Here, partly for predictability and simplicity and partly in order to follow 
common assumption that we live at a typical place in our galaxy, we take a common boost 
factor for all particles in our calculations. The mUED agrees well in the low energy region 
E < 10 GeV, but starts deviating significantly from the recent PAMELA data when energy 
of antiproton becomes higher, if MED diffusion model is adopted. Therefore, mUED seems 
to be disfavored by PAMELA data in the high energy region, and it agrees with PAMELA 
observations only when MIN diffusion model is chosen. However, one can easily see that 
how the situation can be improved in the split-UED case. By enhancing the mass of q± to 
be twice of the li, the prediction of p/p is reduced by more than a factor of two, and will be 
much suppressed if q\ becomes heavier, as shown by the blue and magenta lines in FigJHl We 
also note that if we adopt MIN model, the split-UED case agrees with PAMELA very well 
even for m qi = 2m^ . The lower bound of the mass of q± in split-UED can be estimated by 
the upcoming data of PAMELA in the higher energy region, since a small bump is predicted 
at E ~ 200 GeV, and the height of deviation from background is controlled by the mass of 

C. Gamma-ray 



For diffusive gamma-ray, there are galactic and extragalactic contributions from the an- 
nihilation of LKP dark matter B\. The flux of the gamma-ray from the extragalactic origin 
is estimated as 46] 
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) dJ y 

~dE 



B F E 2 ctf Bi p, 



eg 4:7im 2 Bi H ^l 
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< av >i 
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dy 



d((l + z)E) 



(1 + z)- 3 / 2 e~ z / z " 



1 + rfe( 1 + ^) 



(30) 

where &m and are the density parameters of B\, matter (including both baryons 
and dark matter) and the cosmological constant, respectively; p c is the critical density; H 
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FIG. 8: Antiproton to proton ratio predicted in mUED and split-UED models, together with 
experimental data. In mUED, both MED (red line) and MIN (green line) diffusion models are 
shown. For split-UED, the cases of m qi = 2mj 1 (blue line) and m qi = (magenta line) are 

presented, using MED diffusion model. 

is the Hubble parameter at the present time; z is the redshift, and z up = m^jE — 1 since 
the maximum energy of photon in the rest frame of annihilation is E = m^; the term 
e~ z/2m " with z max (E) ~ 3.3(£7lOGeV)~°- 8 takes the optical depth into account [46|. For 
the numerical results, we use |47| 



Q Bl h 2 = 0.1099, Q M h 2 = 0.1326, Q A = 0.742, p c = 1.0537 x 10" 5 GeV/cm 3 . (31) 

On the other hand, the gamma-ray flux from the annihilation of B\ in the Milky Way 
halo is 



dE 



halo 



E 2 1 diV 7 
~kK2m 2 Bi dE 



p 2 (i)d£ 



l us 



(32) 



where phaio is the density profile of dark matter in the Milky Way, (^j los p 2 {l)ditj is the 
average of the integration along the line of sight (los). 

Since there are still big uncertainties and ambiguities for modeling the Galaxy center, 
we integrate over the whole sky except for the zone of the Galactic plane (i.e. exclude the 
region with the galactic latitudes \b\ < 10°). For the background, we use a power- law form 
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FIG. 9: The diffusive gamma-ray predicted in mUED and split-UED models, together with 
EGRET data. The solid lines are for signals with background, blue line is for mUED; green, 
magenta and red lines are for m qi = 2mi x , m qi = 3mj 1 and m qi — > oo, respectively. Signals of dark 
matter annihilation are shown separately in dashed lies for leptonic final state (red), hadronic final 
states in mUED (blue) and in split-UED with m qi = 2m; 1 (green) and m qi = 3m; 1 (magenta). 



adopted in Ref. [48] 



dE 



5.18 x 10~ 7 E' 
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-i 



(33) 



bg 



where E is in units of GeV. 

The total flux of diffusive gamma-ray is shown in Figl9], and the signals from the hadronic 
and leptonic final states in B\B\ annihilation are presented as well. We shou 



that we do not try to fit EGRET data 



49, 



d emphasize 
50| in the plot, since Fermi 13j is expected 



to have more precise results soon, however, we still show EGRET data for reference. We 
notice that the extragalactic contribution from cosmological distance is very small even we 
enhance it by a factor of 10 6 to estimate the effect of subhalo 46(. The gamma-ray in the 
decay of ir° from hadron cascade decay is softer than that from the decay of r, therefore 
the distribution become flater when more quarks are produced in Bi annihilation. By 
suppressing the fraction of the quark mode in the final state, i.e. comparing mUED (blue 
lines) and split-UED (green, magenta and red lines), the starting point of deviation from 
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the background will shift to high energy 6 and the predicted peak at energy about 200 GeV 
will become smaller. The bump at the high energy is a generic prediction for split-UED 
even all of the higher order KK quarks are decoupled, because gamma-rays still come from 
the r decay, and it can be further examined by the Fermi experiment in the near future. 



IV. COLLIDER SIGNATURE OF SPLIT-UED 



Having light colored particles below a few TeV in split-UED, the LHC can produce lots 
of KK quarks and gluons via QCD interactions. As it is already shown in Sec. [TTJ one of the 
features of our model is that KK quarks are split from the other particles, therefore, will lead 
to collider phenomenology quite different from the mUED models. For example, the 1st KK 
quark, where we take m qi ~ 2m b 1 through calculations, mainly decays to the 1st KK gluon 
gi and a SM quark, generating at least one high p T jet due to the big mass gap between 
q± and g\. Another one is the existance of tree- level KK number non-conserved interactions 
between KK-even gauge bosons and the SM fermions, as a result, the production cross 
sections of 2nd KK gauge bosons can be substantial. It has been shown that the signals of 
a 2nd KK gluon production can be dominant over the SM background in dijet events after 
imposing a certain set of cuts at the LHC [9|. In this section, we will focus on the scenarios 
of 1st KK quarks and gluon, which can only be produced in pairs. 

The production cross sections involving q± at the LHC depend on the bulk mass parameter 
[i. Among the production cross sections of qiq^p, qiq\ and q\q(\ the cr(qiqi) is dominant 
one as the valence u and d partons can contribute to the initial states. These productions 
proceed mainly through the t-channel diagrams with the 1st KK gauge bosons exchanged. 
The most relevant couplings for these productions are qi-q-Vi, where V\ — g±, B\, the 1st 
KK W— boson W\ and the 1st KK Z— boson Z\. The cross section a{uxU\), a{u\di) and 
a{d\di) are plotted in Fig. [10] with varying fi, and we neglect the qiqp productions, since 
they are insignificant. For small fi, the production cross sections increase with increasing 
/i due to the enhancement from the couplings, as shown in Figj3j since the cross section is 
proportional to ({?no) 4 - For large /x, on the other hand, those cross sections decrease very 



6 The preliminary result of Fermi 51 1 shows the data is consistent with the estimated background up to 
E ~ 10 GeV, so by enhancing the mass of qi in split-UED seems to be prcfcrcd, if the preliminary result 
is thoughtful. 
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FIG. 10: The cross sections and the mass of u^i as functions of the bulk mass fi with fixing 
1/R = 620 GeV. Here, u\{d{) includes both uli and uri (dn and dm) contributions. 

quickly with increasing /i (increasing m qi ) as expected. For the giqi and gxgi productions, 
the cross sections, which are not shown, are at the same order as that of q\q\. Here all cross 
sections are calculated using CalcHEP 52j by modifying the mUED model file in Ref. 53]. 

The signature of split-UED qiqi productions followed by q± — > g\q — > B\qX is two 
high pt jets plus soft jets/leptons and large missing momentum. The kinematics and the 
decay branching ratios are very similar to those of the corresponding supersymmetric model 
(SUSY) processes, namely the qq productions followed by q — > gq — > jQqX with squark mass 
m q — m q 1 i gluino mass m~ g = m gi and the lightest neutralino mass m^o = ttib 1 ■ Therefore, 
we mimic split-UED collider signatures of q\ pair productions at the LHC by scaling the total 
cross sections of SUSY signatures with the model point whose mass spectrum is similar to 
that of the split-UED. Furthermore, because collider signatures of q\q\ are the same as that 
of qq, they share the same SM background. Since the SM background for squark and gluino 



productions at the LHC is very well studied 54|], we apply the same cuts in our calculation. 
The model points we study are listed in Table IIV| in which we use ISAJET 55[ for SUSY, 
and we take 1/R and /x to be 620 GeV and 700 GeV, respectively, for split-UED so that 



9i 



2m b 1 - The events are generated by using HERWIG [56|] for SUSY processes, then we 



regard them as the split-UED events, and detector simulations are carried out by AcerDET 
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split-UED 


mass 


SUSY 


mass 


qli 


1347 GeV 


ul Al 


1355, 1358 GeV 


uri 


1322 GeV 


ur 


1304 GeV 


dm 


1318 GeV 


d-R 


1263 GeV 


91 


794 GeV 


9 


799 GeV 


Bi 


621 GeV 


X? 


622 GeV 



TABLE IV: Mass spectrum of split-UED at 1/R = 620 GeV and \t = 700 GeV. The corresponding 
SUSY point generated using ISAJET is also listed. 





after standard cut 


M eff > 1 TeV 


M eff > 1.5 TeV 




0.40 


0.37 


0.21 


1\9\ 


0.30 


0.18 


0.049 


9\9\ 


0.18 


0.04 


0.007 



TABLE V: The acceptances for the split-UED q\q\, q\9i and g\g\ productions after imposing basic 
event selection cuts f)34f) and harder M e ff cuts for the model point in Table IIVI Note that these 
numbers are obtained from SUSY events. 

[57( 1 . Table [V] shows the acceptances of the events after imposing the following basic cuts: 

M eff > 500 GeV, E Tmiss > max(100 GeV, 0.2M eff ), n 100 > 1, n 50 > 4, (34) 

M e ff is the scalar sum of the pr of the first four leading jets; -ETmiss is the transverse missing 
momentum; nioo (^50) is the number of jets with pt > 100 (50) GeV. Forthermore, we also 
use harder cuts on M ef[ with M e ff > 1, 1.5TeV. 

We can see that the acceptance is about 40 % for q 1 pair production after the basic cuts. 
On the other hand, the efficiency to select the g\ pair production is much smaller, especially 
when a large M e g is required. For example, the number of events from the 1st KK gluon 
pair production reduces by factor of 1/30 after imposing the basic cuts with M e s > 1 TeV. 
This is because the probability of having high px jets and large -Ermiss is very low for the g\ 
production events at our point. Thus, we expect g\g\ production is not promising due to the 
very small efficiency found in the simulation. The g\q\ production may be easier to detect 
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FIG. 11: M e g distributions for "gigi" , "<7i<7i" and "qiqi" from left to right obtained by correspond- 
ing SUSY channels gg, qg and qq. We neglect the y-axis here since we focus on the shapes, and we 
see that the distributions shift to higher energy region from left to right. 

compared to the gigi case with the help of the high jet from q\ —>■ qgi, although the signal 
and the background separation would be worse compared with qiqi production. Therefore 
we only consider qiqi productions in the following studies. For completeness, we show the 
M eS distributions for gigi, q\g± and q\q\ in Fig. [TTJ Again these results are obtained from 
SUSY events mimicking split-UED with the same kinematics. 

The total production cross section of qiqi, including uliAli, u Rii dm, is 7.64 pb, and we 
expect 7640 events for 1 fb _1 . According to Table [TTl we expect 7640 x 0.37 = 2830 events 
left under the standard cut with ALs > ITeV. We should also consider lepton veto before 



using the background studied in Ref . 54j . Although it depends on the lepton branching ratio, 
we expect more than a half number of events, i.e. 1400, remain with the lepton veto. The 
number of SM background events with the same cut for lfb -1 is less than 300 according to 



the M e fj distribution shown in Ref. [54|. Therefore the qiqi signal distribution is well above 
the background distribution for M c g > 1000 GeV. 

Given the enough statistics, we now discuss the possibility to determine some of the 



split-UED particle masses. Fig. [12] shows the M?2 |58| distribution of the two highest pt 
jets under the standard cut with M c g > 1000 GeV. The M?2 is usually used to determine 
the masses of unknown particles which are pair-produced and decay identically. In our case, 
we have two gi's being produced, and both decay into g±q. Here the Mt2 is calculated 
from the two highest px jets, j\ and j'2 and a missing transverse momentum defined as 
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FIG. 12: The Mj<2 distribution for q\q\ production events under the standard cut with M e g > 1 
TeV for 1 fb _1 . We can see that the end point of Mt2 distribution is given by m gi — ~ 880 GeV. 
The bin size is 10 GeV/bin. Note that we do not apply any lepton veto cut for this plot. 

PTmiss = ~PTj 1 — Prj 2 and the formula is 



where qr\ and qx2 are two dummy parameters that make up prmiss; and the minimization is 
taken for all possible sets of qn and qri\ Vhkh) * s ^ ne momentum of ji(2); Atrial is the trial 
mass that represents the unknown mass of the daughter particle (gi in our process) and My 
is the transverse mass. The Mt2 depends on the trial mass m tr i a i, and we took m tr i a i = 
in Figfl2l For our simulation, the two highest px jets comes from qi — > giq, and the end 
point of Mt2 should be m qi when the trial particle mass m tria i is taken as the true mass of 
g-y. Another important fact is that when there are no initial state radiations, the formula of 
the Mt2 endpoint for m tr i a i = has a general form given as 



M T2 



mm 

PTmiss=9Tl +?T2 



[max{Af T (4n,Pj 1 , mtriai), M T (q T2 ,p j2 , m tria i)}] , 



(35) 



Mr 




(36) 
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where ttia is the mass of mother particle which is pair produced and mx is the mass of the 
unknown daughter particle, and in our case rriA = rn qi and mx = m gi . We can see that 
the Mt2 distribution in Fig. [T2] has a clear end point which is consistent with the predicted 
value M^ d — 880 GeV for the model point in Table HVl . Although we are unable to measure 
precise masses of q\ and g±, we still have a useful information on their combination. 

Finally, let us comment on the comparison with mSUGRA. If we assume the signature 
comes from squark production, with a typical mass spectrum rriq ~ 1 TeV and gives the 
similar M?2 end point, the cross section is much smaller than that of the split-UED point. 
Thus, SUSY and split-UED should be distinguished based on the event rates, although the 
kinematical nature of the signal is similar. We can get a rough understanding on the enhance- 
ment of the cross section in split-UED by comparing the helicity structure of squarks / gluinos 
and corresponding KK quarks/gluons 



7 See also 59 1 for similar phenomenon in Little Higgs model. 
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V. CONCLUSION 

The split-UED model proposed in Ref. j^l is studied in detail. From the basic set-up 
of the model, we calculate all the masses (including one-loop radiative corrections to the 
masses) and relevant couplings for the first KK particles. Characteristic split spectrum of 
KK quark states is realized and enables us to control the hadronic branching fraction of the 
LKP annihilation by one parameter, bulk mass scale (/i). In the non-relativistiv limit, such 
a branching fraction decreases very fast by increasing fi. 

In order to explain the anomalies from the cosmic ray, we fixed the size of extra dimension 
(1/R — 620 GeV) from the peak position in flux of electrons plus positrons observed by ATIC 
and PPB-BETS, with boost factor B = 200. Then all data from PAMELA, ATIC, and PPB- 
BETS can simultaneously be fitted very well. It is important to notice that such a dark 
matter mass in split-UED generically can explain the right amount of relic density of dark 
matter (Q ~ 0.23). To avoid all the strong constraints known up to date, we take m qi >2/R 
so that the hadronic branching fraction of LKP annihilation is significantly smaller than the 
leptonic one {pqqjou < 10%) and the cosmic antiproton flux is suppressed more than the 
safe rate. The flux of cosmic gamma-rays from pion decay is also highly suppressed and 
hardly detected in low energy region (E^ < 20 GeV). However future coverage of high 
energy domain in O(100) GeV will make it possible to detect gamma-rays originated by tau 
lepton which is still sizable. 

At the LHC, the 1st KK colored particles (KK quarks and gluons) in split-UED are 
copiously produced and further decay into hard jets, soft jets/leptons and dark matter 
(missing energy). Because of the large mass splitting between the 1st KK quark qi and 1st 
KK gluon gi, we can separate out the production of gigi, g\q\ and q\q\ by using the M e g 
cuts. We focus on the q\q\ production, in particular, which has two hard jets so that we can 
distinguish the signals from the SM background. The split-UED q\ pair signal is simulated 
by scaling the SUSY signatures with the same mass spectrum. After the appropriate cuts, 
our signal is well above the SM background. We calculate the Mt2 distribution of the 
qi pair signals, and its end point reflects the information for combination of q\ and gi 
masses. Although we are unable to determine the mass of individual particle in split-UED, 
its predictions of the signatures of q\qi production at the LHC, e.g. M e g and M^ 2 , can be 
examined whether they are consistent with the cosmic-ray signatures in the near future, and 
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vice versa. 
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Appendix A: KK decomposition in split-UED 



For all the gauge bosons, scalars, and leptons, their KK decompositions are the same as 
those in mUED. For bosons and fermions that are in the same chirality of the zero mode, 
which has (+, +) boundary condition at [—L,L], the KK decomposition is: 

rTix 



=(f) (x) h — j= 

2L \[L 



sin 



2L 
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(x) + cos 
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2L 



y 



(37) 



While for fermions that are in the opposite chirality of the zero mode, which has the (— , — ) 
boundary condition at [—L,L], the KK decomposition is: 



rr 
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n T ,n 



n it 



2L 



i n - (x) + sin 



2L 



<p n +{x) 



(38) 



The label n~ or n + here stands for the n-th KK modes with the even/odd KK parity 8 . 

For the quark, we consider the case in which a SM quark is embeded into \P + = P + ^> 
component of a 5D Dirac fermion where P + = (l+7s)/2 is the positive chirality projection 
operator. The other case could be considered quite similarly by re pla cing \x by A 
convenient way of expressing the KK decomposition is in the form [9!. Il^: 



V+(x,y) = 9n+(\y\)Xn+(x) + e(y)g n -(\y\)x n -(x), 



= Yl <y)M\v\)Mx) + fn-(\y\)lpn-{x). 



(39) 



n T ,n 



The 5D profiles satisfy the following coupled, first-order equations of motion 

d y g n + H9n - rn n f n = 0, 

9yfn ~ fJ>fn + m n9n = 0, 



(40) 



In our convention, the KK number n = 2n — 1 and n + = 2n + for the {n , n + } in Ref. 14] 
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with each 5D profile g n +, g n ~, f n +, f n - satisfying the (+, +), (-, +), (-, -), (+, -) boundary 
condition at y — and L, respectively. Once the solution for y C [0, L] is obtained, the 
solution to the whole space is determined from Eq. (I3"5j) thanks to the symmetry. 

For the zero modes, m n = 0, the equations are separable and the solution is given as 



9o+(y) = A exp( / \{$(s))ds) = A exp(iiy), 



(41) 



with the normalization factor A = \J yu/e^ 2 ^ -1 '. There is no zero mode company for g + 
but each of the KK modes has its couple and fill the spinor states of the Dirac spinor. The 
solution for the even KK modes is 

k r , 



9n 
fn 



cos(k n +y) 



L lm % 



m r 



sm(k n +y) 



(42) 



where the KK mass m n + = \/ [J- 2 + k^ + and k n + = nir/L (n G Z). 

To avoid the very light 1st KK quark, we choose [i > and the zero mode is quasi-localized 
at the boundary y = ±L. For odd KK modes, the solution is 



9n~ 
fn~ 



-sin(k n -y), 

\ \^— cos(k n -y) + sin(k n -y) 
L Lm„- m„- 



(43) 



where the KK mass m n = ^//x 2 + k\ and k n is the n-th solution of the equation 



k r , 



-/itan(fc n -L) . 



(44) 



When [i increases from to +oo, k n - increases from (n — 1/2)tt/L to mi/ L. In this case, in 
the limit of fi — > +oo, all KK modes could be decoupled. 
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